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Abstract. Recently we have studied, using numerical simulations, a limited valency 
model, i.e. an attractive square well model with a constraint on the maximum number 
of bonded neighbors. Studying a large region of temperatures T and packing fractions 
cj>, we have estimated the location of the liquid-gas phase separation spinodal and the 
loci of dynamic arrest, where the system is trapped in a disordered non-ergodic state. 
Two distinct arrest lines for the system are present in the system: a (repulsive) glass 
line at high packing fraction, and a gel line at low <p and T. The former is essentially 
vertical (^-controlled), while the latter is rather horizontal (T-controlled) in the (4>—T) 
plane. We here complement the molecular dynamics results with mode coupling theory 
calculations, using the numerical structure factors as input. We find that the theory 
predicts a repulsive glass line — in satisfactory agreement with the simulation results 
— and an attractive glass line which appears to be unrelated to the gel line. 



1. Introduction 

In recent years, investigation of structural arrest in colloidal systems has witnessed 
a renewed and large interest in the scientific community since new phenomena have 
been identified, in particular when the particles interact via an attractive potential 
of range short enough compared to the size of the colloid PQ. Theoretical El El, 
numerical 00 Ej and experimental 00110111111121113 studies have shown the existence 
of two different mechanisms responsible for the slowing down characteristic of structural 
arrest. At high temperature T and high packing fractions <ft of the dispersed phase caging 
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effects prevail and produce the typical repulsive glass behavior. At low T, and slightly 
smaller <ft, another mechanism sets in, due to the stickiness of the particles, and generates 
a so-called attractive glass [Hj. It was hypothesized at the beginning that the latter 
mechanism could be one of the routes to the formation of a gel at low densities j^j. 
Many investigations of these phenomena have been made fSl ED- but most of them 
were faced with difficulties related to the existence of a two-phase region at low <p. In 
fact it was unambiguously shown, for the case of short range attractive potentials, that 
the arrest line intersects the binodal line at its high volume fraction side [T7] ITH] . The 
crossing may be avoided if a lowering and shrinking of the two-phase coexistence region 
is achieved. 

One possibility to limit or suppress phase separation is to consider the effects 
of a long-range repulsion complementing the short-range attraction. This type 
of potential, able to mimic effects of screened electrostatic interactions, properly 
describes interactions in charged colloidal suspensions. This route has been actively 
pursued very recently in experiments (THl EH1 EI] simulations [221 EH 1211 123 EE] and 
theory [23 EE1 El EH • 

Another possibility to limit or suppress phase separation, which does not invoke 
the presence of repulsive long range repulsion is offered by saturation of bonding, i.e. 
by a limit on the maximum number of bonded interactions. To prove this mechanism, 
we have devised a model where an ad-hoc constraint is adopted, in addition to a square 
well (SW) interaction. We studied numerically a saturated SW model first introduced 
by Speedy and Debenedetti [21 E21 ESj - also named limited valency model. In this 
model, the square well attraction between two particles is constrained to a maximum 
number of bonds n max that particles can form. By reducing n max down to low values 
such as 3,4,5, we found that the location of phase separation is progressively shrunk 
both in <p and T [33]. Hence, reduction in the number of bonds may avoid the crossing 
between the dynamic arrest and the binodal lines, or to move the crossing to much 
lower and T. Simulating this model, we showed that indeed the avoidance of phase 
separation allows the emergence of arrested states at low 0, which have quite different 
physical features than glasses, both of attractive and repulsive type. Completing the 
study up to very large densities in [HE], we detected the presence of two distinct arrest 
lines: a gel line, practically flat in T , governed by Arrhenius dynamics, and a glass line, 
practically constant in 0, corresponding to the standard hard sphere glass transition. All 
along the gel line, the dynamic features are not reducible to those of attractive glasses, 
suggesting that the two states in attractive systems do not necessarily correspond to 
the same phenomenon, and a simple extension of the attractive glass line may not be 
always appropriate. 

The purpose of the present work is to use the MCT in order to check the conclusions 
obtained through numerical simulations, to clarify in particular the possibility of 
describing the low <fi equilibrium gel formation. To this aim, we compare here our 
numerical results of Refs. [331 EH] with predictions of Mode Coupling Theory (MCT) for 
the same model, using as input of the theory the "exact" structure factors for the n max 
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model, calculated numerically. 

2. Overview of the simulation results 

The limited valency model described above has been recently studied in detail and we 
summarize here the main results of the numerical investigation |3*4^ I3*H] . Molecular 
dynamics was performed in a square-well attractive system of N = 10 4 particles of unit 
mass and hard core a = 1 with a range of the well A such that (a + A)/er = 1.03 
and depth Uq. Temperature is measured in units of uq, with the Boltzmann constant 
ks = 1- The maximum number of bonds was limited to n max = 3,4,5 as compared to 
the unconstrained SW, where n max = 12. Indeed, for n max > 6, no significant difference 
in the phase behaviour was observed with respect to the standard square well model. 
The main findings are reported in the following list summarizing our results. 

(i) The coexistence region in the (0, T) plane reduces in size as n max decreases, the 
values of above which the spinodal disappears are ~ 0.20 and ~ 0.30 for 
n max = 3 and 4 respectively. The loci of percolation also shift to lower T values, 
and cross the two-phase region on the low side, as shown in Fig. ^ 

(ii) The bonds formed due to the attraction show a lifetime which follows an Arrhenius 
behavior in T . 

(iii) The mean squared displacement develops a plateau at low T which defines a 
localization length of the order of the particle size, much larger both than the 
length typical of a repulsive glass due to caging (of the order of 0.1a) and than 
that of an attractive glass (of the order of the well width). The diffusivity shows 
a power-law behaviour in along constant temperature paths, while it displays an 
Arrhenius dependence on T for constant volume fraction. 

(iv) The normalized intermediate scattering functions show a corresponding plateau 
at low T only for values of the momentum transfer q smaller than the value 
corresponding to the first peak of the structure factor. This feature is rather 
different from the standard behavior close to a glass transition, and is similar to 
the behavior of a chemical gel 

(v) The height of the plateau for the normalized intermediate scattering functions, i.e. 
the non-ergodicity factor f q , has a different shape in q for the gel than both its 
attractive and repulsive glass counterparts. 

Based on these facts we concluded that the system shows a glass line at high and a 
gel line at low T and low 0. More important, the dynamical behavior observed close to 
the gel line does not seem to allow the identification of the gel with the attractive MCT 
glass. 



Mode- coupling theory predictions for a limited valency attractive square-well model 4 

T 

0.4 



0.2 



0.0 

0.0 0.1 0.2 , 0.3 

Figure 1. Spinodal and percolation lines (dashed and solid lines respectively) of the 
limited valency model for n max = 3 (circles), n max — 4 (squares), n max = 5 (triangles) 
and n max = 12 (no symbols). 

3. Mode Coupling Theory 

We solve the MCT equations to locate the ideal glass line(s) in the (0 — T) plane. The 
MCT, starting only from the structural information contained in the static structure 
factor S(q), provides indication of the onset of non-ergodic behavior [37 . Although 
the theory is not strictly based on the hypothesis of pair additivity for the interaction 
potential, its main (uncontrolled) approximation is the factorization of higher order 
correlation functions into products of pair correlation functions [37J[2H|- Since the n max 
model incorporates many-body terms in the Hamiltonian |34[ 13*3] . we cannot expect the 
theory to work at its best. Perhaps higher order correlations should also be considered, as 
in the case of silica, a tetrahedral network-forming liquid, where the triplet correlation 
function C3 provides a relevant contribution to the MCT kernel "2~IJ- Here, we limit 
ourselves to the standard version of MCT calculation, because the random nature of the 
^max bonds along the particle surface does not produce any angular constraint, retaining 
the full sphericity for the model. However, it would be interesting to check in a future 
work whether the inclusion of the corresponding three-point correlators calculated within 
MD simulations produce any significant difference in the MCT results. 

MCT is able to predict the full time evolution for the density- density autocorrelation 
function $(g, t) as a function of the momentum transfer q and t, through coupled 
non-linear integro-differential equations. For a given interaction potential, through 
the knowledge of S(q), the memory kernel entering in the nonlinear term of the MCT 
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equations can be evaluated and the equation solved for various values of q. In particular, 
the non-ergodicity transition leading to structural arrest is obtained by performing the 
limit t — > oo in the equations. Defining the non-ergodicity factor as the long-time limit 
of $ correlator lim^oo t) = f(q), f(q) is found to be the solution of [37], 

where the memory kernel m(q) is quadratic in the correlator itself, 

1 f d 3 k 

m(q) = -j (2^5 V(q, k)/(*)/(|q - k|). (2) 

The vertex functions V, the coupling constants of the theory, are 

V(q, k) = ^ [q • (q - k) c(|q - k|) + q • k c(k)} 2 S(q)S(k)S(\q - k|) (3) 

and depend only on the Fourier transform of the direct correlation function c(q), or 
equivalently on S(q), and on the number density p. In the A2 bifurcation scenario 
of MCT |2Zj the solutions of Eq. jump from zero to a finite value at the ideal 
glass transition. The locus of the fluid-glass transition can be calculated varying the 
control parameters of the system, <ft and T. For a square well model, there is an 
additional control parameter, that is the range of attraction A. In this case, higher 
order bifurcations of the solutions arise when A < A* ~ 0.041 [4 j. In this case, two 
distinct glassy solutions appear, a repulsive and an attractive glass respectively, with 
different non-ergodicity parameters and mechanical properties |4()j . 

We calculate the MCT ideal glass transition line for n max = 3 using as input the 
numerical S(q), "exact" within numerical precision. In the evaluation of the MCT kernel 
of Eqn. (j2J, it is crucial to integrate over all q contributing to the memory function. 
For short-range attractive potentials, it is important to integrate up to very large q 
values, since the information of the potential shape is coded into the large q region. 
From a numerical point of view, it is convenient to evaluate S(q) at large q by Fourier 
transforming the pair distribution function g(r) and at small q, by direct evaluation of 
S(q) in g-space. Indeed, bonding in g(r) is reflected by a very large and constant signal 
between a and a + A, which fully accounts for the large q behavior of S(q). For the 
same reason, the large-g 'beats' |Hj in S(q) have the same shape as the standard SW 
and carry information of the short-range potential. The asymptotic form of S(q) for 
large-g, both for the SW and the n max model, is 

SasympM) = 1 - ^ {sin (qa) - qo cos (qa) 

+ (e^ - 1) [q(a + A) cos (q(a + A)) 

— sin (q(a + A)) + sin (qa) — qa cos (<?cr)]}. (4) 

The amplitude A (which depends on T and <ft) is different for the SW and n max models, 
being smaller in the n max case, due to the reduced number of bonded neighbors. 
Combining information from S(q) and g(r), an accurate description of S(q) over the 
entire relevant q range is obtained, as shown in Fig. El In the SW case, the large-g 
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Figure 2. Inset: Static structure factor for <f> = 0.20 and T = 0.2 used as input 
in MCT calculated directly from simulations (points), and by Fourier transform of 
the g(r) (line). Main: Enlargement of the large-g behavior displaying 'beats'. For 
comparison we also report the corresponding asymptotic SW (differing only by a factor 
A in amplitude) and HS (out of phase) results. 

oscillations resulting from the narrow width of the square well potential are largely 
responsible for the MCT attractive glass transition [4j. We solve the MCT equations 
in g-space in the window from to 600a _1 , with a mesh of about 0.3<t -1 , for a total of 
2000 q- vectors. We bracket the MCT arrest line by locating two adjacent state points 
along isochores where a liquid and glassy solutions are respectively found. The resulting 
MCT predictions are reported in Fig. El 

At high densities, MCT results are very similar to those found for the simple 
SW jU HO] for the A = 0.03 case. Even in the n max case, the MCT equations predict 
two distinct glass lines, respectively an attractive and a repulsive one, and a glass-glass 
transition ending in an A 3 singularity [42J. A comparison between the MCT predictions 
for the SW and for the n max model is reported in Fig. El The glass lines for both models 
converge to the hard-sphere result at high T and to the same locus for <fi — ► 0. This is 
consistent with the expectation that at high arrest is driven by packing, while at low 
the constraint on the maximum number of bonded neighbors becomes irrelevant and 
hence the two models tend to become similar. Differences with respect to the SW case 
are observed at intermediate 4>. 

In the SW case the ideal MCT attractive glass line is rather flat, monotonically 
decreasing in T, almost merging into the spinodal on the left side of the critical point |40j. 
The n max model shows instead a non-monotonicity of the T-dependence of the attractive 
line with decreasing 0, causing the presence of a minimum in T for <fi « 0.4 and T « 0.24. 
The significant suppression of the attractive glass line in the n max model at intermediate 
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Figure 3. Comparison of MCT predictions and spinodal lines for the SW case and 
for the n max model with n max = 3. The inset shows an enlargement of the high (f> 
region, to visualize the attractive and repulsive lines and the glass-glass transition. 

arises from the decrease in the number of bonded nearest neighbors as compared to 
the SW case. The attractive glass line has a maximum around ~ 0.10 and T « 0.3, 
before turning down following the spinodal as in the SW case, on the left of the critical 
point. It is possible that the shape for < 0.1 in the MCT line is an echo of the 
underlying increase of the compressibility at low due to the close-by spinodal, since 
at low the MCT line roughly follows in shape the loci of constant 5(0) (see Fig. |UJ). 

In the SW, the attractive and the repulsive glass lines differ essentially in the q- 
dependence of the non ergodicity parameter f q . The theory also allows us to calculate 
f q and its critical value at the glass transition. As compared to the repulsive glass, 
the attractive glass is characterized by much larger f q values, extending to much larger 
q. These features are also displayed in the theoretical calculations for the n m3X model. 
Along the attractive glass line, for all values, from about = 0.05 up to = 0.54, 
the critical non-ergodicity parameter, i.e. f q at the MCT transition, does not change 
significantly. Fig. |U shows the full g-vector dependence for both types of glasses. We 
have also verified that in the attractive glass phase f q is significantly dependent on T. 

For comparison, we plot in Fig. |3] the non-ergodicity parameters estimated from 
the simulations, via a stretched exponential fit of the correlators jHH], along the iso- 
diffusivity line DT~ l l 2 = 0.0005. A transition from a gel to a glass is evident. However, 
the shape of f q never resembles that of an attractive glass. The gel is characterized by 
non-ergodic features only at large length-scales, due to the mobility of the network in 
the available free space, while an attractive glass is strictly confined within the bond 
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Figure 4. Critical non-ergodicity parameters f q calculated within MCT. Note that 
at low <j) along the attractive glass line, f q are indistinguishable. 




Figure 5. Critical non-ergodicity parameter f q calculated from the simulations (i.e. 
along the iso- diffusivity line DT^ 1 / 2 = 0.0005), via a stretched exponential fit of the <f> 
correlators. Note that f q never resembles that of an attractive glass. Lines are guide 
to the eye. 
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length. The main difference between theoretical and simulation results is that, although 
also in the gel bonds are on average permanent on the time scale of the simulations 
at low T [HE], the limited number of neighbours accounts for residual motions of the 
particles. Hence, particles are confined by the attractive well width only relative to each 
other, but still can freely explore up to their diameter length without ever breaking the 
network. In the theoretical calculations, however, the presence of the short-range bonds, 
manifested in the large-g tail of S(q), is responsible for the non-ergodic transition. 

Already at the SW level, the theory strongly overestimates the tendency to form 
an attractive glass (4*3*1 I44j . Compared with numerical simulation (2*3, the theoretical 
attractive glass line has to be shifted approximatively by a factor of 3 to 4 in T, while 
the repulsive glass line needs only a 10 per cent adjustment in 0. The correction of the 
location of the ideal MCT lines provokes a major effect: the attractive glass line does 
not lie above the critical point but meets the spinodal at low T on the right side |17j . 
If the overestimate of the attractive glass line is properly taken into account, then one 
has to conclude that it is not possible to form arrested states at low without the 
intervention of a phase separation. Recent simulation studies confirm that this is the 
case independently from the width of the attraction range P-8 . It is also necessary to 
recall that the ideal MCT glass lines, especially when energetic caging is dominant, have 
to be interpreted as cross-over lines, from a power-law to an activated dependence of 
dynamic properties (20! ■ We can summarize the dynamical arrest behavior in Fig. 101 
One locus of arrest is found at high 0, rather vertical and corresponding to the hard- 
sphere glass transition. This locus is quite well described by MCT. Very different is the 
situation concerning the low T slowing down. The isodiffusivity lines suggest a rather 
flat arrest line. Two different loci could be associated with arrest at low T. One defined 
by the T c of the power-law fits of the diffusivity and one at T = associated with the 
vanishing of D according to the Arrhenius law [3"*"]. It would be tempting to associate 
the T c -line with the attractive glass line predicted by MCT and interpret the wide region 
between the two lines as a region of activated bond-breaking processes [43J. 

Results for the present model provide evidence that a gel line can not necessarily 
be considered an extension of the attractive glass line. Here, a gel never transforms into 
an attractive glass, but only into a repulsive glass with increasing 0. The competition 
between the two arrested states seems also to produce anomalies in the dynamics of 
the same kind as those found in the presence of MCT higher order singularities |35). 
In particular, these anomalies correspond to a logarithmic relaxation for the density 
correlators in the liquid region close to the singularity, accompanied by a sub-diffusive 
behaviour for the mean squared displacement (45 1 . 

4. Conclusions 

The aim of this work was to compare results for the dynamical arrest in the n max 
model for n max = 3 with MCT predictions. Using the numeric S(q), "exact" within the 
precision of the statistical averages, we have solved the MCT equations and evaluated 



Mode- coupling theory predictions for a limited valency attractive square-well model 10 
0.5 

0.4 

0.3 

T 

0.2 
0.1 
0.0 

0.0 0.1 0.2 0.3 ^ 0.4 0.5 0.6 

Figure 6. Summary of the thermodynamic and kinetic phase diagram for n max = 3, 
including spinodal (dashed lines with filled circles), iso-diffusivity loci where D/y/T = 
0.005,0.0005 (lines with triangles and squares), iso-S(O) locus where S(0) = 1.5, MCT 
lines (attractive and repulsive), extrapolated glass, labeled as <p c , and gel, labeled as 
D = 0, lines respectively from power law and Arrhenius fits (T c = 0) 35 . 

the glass lines. The theoretical results differ only slightly from the SW case and do not 
provide any indication of a gel line. According to the theory an attractive glass line 
(i.e. with the typical features of the attractive glass as localization length ~ A 2 and 
very wide f q ) should be present. From a theoretical point of view, such a line arises 
from the large q oscillations in S(q), i.e. from the g-space signature of the short-range 
bonding. Hence, it has the same origin as in the SW. In the simulation of the n max 
model the bond localization is not observed, neither in the MSD nor in the width of 
f q . We believe this is due to the fact that, although bonding is present, particles are 
confined by the potential well only relative to each other. Indeed, at low <f), oscillations 
of parts of the connected network are possible, preventing the observation of the MCT 
mechanism for arrest. At small 0, even at extremely low T, when the bond lifetime 
becomes comparable to the simulation time and bonds between particles are essentially 
permanent, it was shown in that the plateau of the MSD remains of the order of the 
particle size. Of course, the inclusion in the MCT calculation of higher order correlation 
functions could be important in the study of a many-body interaction and should be 
considered. However, due to the sphericity of the model, our opinion is that this should 
not make a significant change, i.e. the prediction for the attractive glass transition 
should be robust in this respect. 

In summary, the present study provides a clear indication that in the present model, 
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where the liquid-gas phase separation can be avoided and arrest at low <fi can be explored 
in equilibrium conditions, the observed arrested state is substantially different from the 
the low-0 extension of the attractive glass. In the present model, in which bonding 
interactions are not strongly directionally constrained, large amplitude motions are 
possible within the fully bonded network. These modes significantly affect the shape of 
the correlation functions and make it impossible to observe the short-range localization 
characteristic of the attractive glass. 

The comparison with the solution of the MCT equations for the n max model (a 
possibility offered by the spherical symmetry of the interaction potential) confirms that 
MCT significantly overestimates the role of the bonding, predicting an attractive glass 
even in the present case. The present results strongly suggest that the attractive glass 
is an arrested state of matter which can be observed in short-range attractive potentials 
only at relatively high 0, being limited by the spinodal curve. When the inter-particle 
potential favors a limited valency, arrest at low <fi in the absence of phase separation 
becomes possible but with a mechanism based on the connectivity properties of the 
particle network. The dynamic features of this slowing down, at least in this model of 
geometrically uncorrected bonds, is clearly different from what would be the extension 
of the (attractive) glass line. 
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